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Abstract 

Various shadowing aspects of nuclear parton distributions are discussed in a parton- 
model framework. We explain existing x dependent data of F^/F®, then the model 
is applied to a flavor asymmetry u — d in a nucleus, to shadowing in valence-quark 
distributions, and to nuclear dependence of Q 2 evolution. First, we find that a finite 
nuclear u — d distribution could be possible due to nuclear interactions even in the flavor 
symmetric case in the nucleon (u — d—0). Second, it is indicated that valence-quark 
shadowing could be used for discriminating among various shadowing models. Third, 
we find that nuclear dependence of Q 2 evolution 9[F 2 5n /F 2 c ']/9[ln Q 2 \, which was found 
by NMC, is essentially understood by modification of parton distributions in nuclei. 
However, higher-twist effects in nuclear interactions could be tested by studying the 
details of the nuclear Q 2 evolution. 

1. Introduction: structure function F 2 in nuclei 

Nuclear modification of the structure function F 2 has been investigated extensively. 
In recent years, it became possible to measure the ratio F 2 A /F 2 D in the "shadowing" 
region, namely at small x. Accurate experimental data make it possible to test the 
details of theoretical shadowing models. Among various models for describing F%, we 
study a hybrid parton model with Q 2 -rescaling and parton-recombination mechanisms 
[1]. The rescaling model was originally proposed by considering possible confinement- 
radius changes in nuclei. Although there is still such a possibility of the nucleon- 
size modification, the major part of the EMC effect at medium x could be explained 
by the nuclear binding. In order to understand the seemingly different ideas, Close, 
Roberts, and Ross claim that it is possible to relate two pictures by using factorization- 
scale independence [2]. In calculating a nuclear structure function in the operator 
product expansion, we separate it into two pieces depending on short and long distance 
physics. This separation scale, which is called the factorization scale, is an arbitrary 
constant and a physical observable should not depend on it. The various interpretations 
correspond to different choices of the factorization scale. Therefore, we may view the 
rescaling explanation as an alternative way of explaining the EMC effect to the binding 
explanation rather than as a different one. 

In the small x region, we use the parton-recombination mechanism for describ- 
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ing the shadowing. The longitudinal localization size of a parton exceeds the average 
nucleon separation in a nucleus at x < 0.1. It means that partons in different nucle- 
ons could interact in the nucleus, and the interaction is called parton recombination. 
Its contribution is a higher- twist effect which is proportional to a s /Q 2 . In order to 
explain x dependence of F 2 Ca /F 2 D measured by NMC, it is necessary to take rather 
small Q 2 (Ql =0.8 GeV 2 ) if the KMRS-B type distributions are chosen as input ones. 
The rescaling and recombination mechanisms are used for calculating nuclear parton 
distributions at Qq. Then, they are evolved to those at larger Q 2 . 

We show theoretical results for x dependence of the ratio F 2 Ca /F 2 D [1] together with 
experimental data. In Fig. 1, the solid curve shows theoretical results at Q 2 =5 GeV 2 . 
They are compared with experimental data by NMC and SLAC [3]. We obtain good 
agreement with the data if Qq is adjusted to a small value (0.8 GeV 2 ). The model 
explains the modification at small x (< 0.05) and at large x (> 0.8) in terms of the 
recombinations, the depletion at medium x by the rescaling, and the antishadowing at 
x ~ 0.2 by competition between two mechanisms. We find that our parton model can 
explain the NMC data F 2 A /F 2 D , so that the model is applied to other interesting topics 
in the following sections. 

2. Flavor asymmetric distribution u — d in nuclei 

NMC suggested that the Gottfried sum rule should be violated in 1991. Since then, 
there have been efforts to investigate mechanisms of creating a flavor asymmetric dis- 
tribution u — d in the nucleon. In order to test the NMC finding, Drell-Yan experiments 
for the proton and the deuteron are in progress at Fermilab. On the other hand, there 
exist Drell-Yan data for various nuclear targets, so that some people use, for example, 
tungsten data for investigating the flavor asymmetry. However, we have to be care- 
ful in comparing the NMC results with the tungsten data because of possible nuclear 
medium effects. In order to find whether such comparison makes sense, we estimate 
a nuclear modification effect. It in turn could be found experimentally by analyzing 
accurate Drell-Yan data in the near future. 

We investigate the u—d distribution in the tungsten nucleus [4]. If isospin symmetry 
could be applied to parton distributions in the proton and the neutron, the distribution 
per nucleon becomes x[u(x) — cI(x)]a = —ex[u(x) — d(x)] pro ton without considering 
nuclear modification. It is just the summation of proton and neutron contributions. 
The neutron-excess parameter e is defined by e = (N — Z)/ (N + Z), and it is 0.196 for 
the tungsten ^Wno- According to the above equation, the flavor distribution should 
be symmetric ([u — d]w = 0) if it is symmetric in the nucleon. However, it is not 
the case in the recombination model. In a neutron-excess nucleus (e > 0) such as the 
tungsten, more d quarks are lost than u quarks are in the parton recombination process 
qq — > G in Fig. 2 because of the d quark excess over u in the nucleus. In the case 
of u — d=0 in the nucleon, only qq — > G processes contribute. If it is not symmetric 
in the nucleon [(u - <i)jv ^ 0], qG — > q type contributions become dominant. The 
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flavor asymmetry is given by x[Au(x) — Ad(x)]A = —{w nn — w pp )x[Au(x) — Ad(x)] pp , 
where w nn and neutron-neutron and proton-proton combination probabilities 

(w nn — w pp = e). [u(x) — d(x)] pp is the asymmetry produced in the proton-proton 
combination. Detailed equations of the recombination effects are given in Ref . [4] . 

We evaluate the recombination contributions with the input parton distributions 
MRS-DO (1993) at Q 2 =4 GeV 2 . In the (u -d) N = case, A = is taken in the MRS- 
DO distribution. Obtained results are shown in Fig. 3, where the solid (dashed) curve 
shows the x[Au — Ad] A distribution of the tungsten nucleus with the flavor symmetric 
(asymmetric) sea in the nucleon. In the {u — d) n = case, the positive contribution 
at small x can be understood by the processes in Fig. 3. In the (u — d) N ^ case, 
the q(x)G — > q process is the dominant one kinematically at small x. Its contribution 
to q(x) is negative because q with momentum fraction x is lost in the recombination 
process. However, more d(x) is lost than u(x) is because of the d excess over u in 
the proton ([Aw — Ad] pp > 0). Therefore, the overall contribution becomes negative in 
[Au(x) — Ad(x)]A due to the neutron excess as shown in Fig. 3. In the medium x region, 
the qG — > q(x) process becomes kinematically favorable. Because it produces q with 
momentum fraction x, its contribution becomes opposite to the one at small x. The 
above results are obtained at Q 2 =4 GeV 2 . Considering the factor of two coming from 
the Q 2 dependence, we find that the nuclear modification is of the order of 2%-10% 
compared with the asymmetry suggested by the MRS-DO distribution. 

[Summary] We find that finite flavor asymmetric distributions are possible in 
nuclei {[u{x) — d(x)]A ^ 0) even in the symmetric case in the nucleon ([u(x) — d(x)]N = 
0). According to the recombination model, the nuclear effects on [u(x) — d(x)]A are of 
the order of 2%-10% compared with the one estimated by the NMC flavor asymmetry 
in the nucleon. Because the Drell-Yan experiments on proton, deuteron, and nuclear 
targets will be completed at Fermilab in the near future, it is possible in principle to 
study the nuclear modification experimentally. 

3. Valence-quark shadowing 

It is shown in section 1 that the parton-recombination model explains the F 2 shad- 
owing fairly well if small Qq is taken. However, it is not the only shadowing explanation. 
There exist other models which produce similar results on the F 2 shadowing. We may 
discriminate among various models by predicting unobserved quantities. Gluon shad- 
owing [5] should be an interesting one for future theoretical and experimental studies. 
As another possible quantity for testing the models, we discuss valence-quark shadow- 
ing in this section [6]. We show that it depends much on theoretical ideas by using 
two different shadowing models. The first one is the recombination model with Q 2 
rescaling in section 1. The second is an aligned-jet model [7], which is an extension of 
the vector-meson-dominance (VMD) model. 

The first model prediction is the following. The rescaling produces depletion in the 
ratio i?3 = Va/Vd at medium x as it explains the original EMC effect. Because the 
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rescaling satisfies the baryon-number conservation, the ratio R3 becomes larger than 
unity at small x. Parton-recombination contributions are rather contrary to those of 
the rescaling. The recombinations decrease the ratio at small x and increase it at 
medium and large x. The overall contributions are shown in Fig. 4 by a solid curve 
(model 1). Antishadowing is obtained in this model instead of the shadowing in F 2 . 
This difference is caused by the fact that the rescaling is used only for the valence-quark 
distribution and not for sea-quark and gluon distributions. 

In the aligned-jet model, the virtual photon transforms into a qq pair, which then 
interacts with the target. However, the only qq pair aligned in the direction of 7 inter- 
acts in a similar way to the vector-meson interaction with the target. The propagation 
length of the hadronic (qq) fluctuation exceeds the average nucleon separation in a nu- 
cleus at small x, so that the shadowing phenomena occur due to multiple scatterings. 
In this model, vector-meson-like qq pairs interact with sea quarks and valence quarks in 
the same manner. Therefore, the valence-quark shadowing is very similar to the F 2 one 
as shown in Fig. 4 by a solid curve (model 2) [7]. The shadowing is calculated by the 
aligned-jet model at small x, then the curve is extrapolated into the medium x region 
simply by using the baryon-number conservation. It is interesting to find completely 
different shadowing results in models 1 and 2. 

Next, we investigate whether both results are allowed by existing experimental 
data. There are neutrino data for the deuteron and for nuclei; however, the accuracy 
is not good enough to test the nuclear modification particularly in the small x region. 
Instead, we use existing F 2 jF 2 data and the baryon-number conservation. We assume 
that Va/Vd is the same as F 2 A /F 2 D in the region x > 0.3 because F 2 is dominated by 
valence-quark distributions. The SLAC F 2 Ca / ' F 2 data [3] at x > 0.3 are fitted by a 
smooth curve in Fig. 4, then it is extrapolated into the small x region by considering 
the baryon-number conservation. Because there is no guideline for the extrapolation, 
a straight dashed-line is simply drawn in Fig. 4. As another possibility, the curve is 
smoothly extrapolated by allowing about 6% antishadowing at x = 0.1 — 0.2 (dotted 
curve). The hatched area between the line and the curve is roughly the region, which 
is allowed by the current experimental data. The figure indicates that both results are 
allowed by the experimental data at this stage. 

[Summary] The valence-quark shadowing is investigated in two different models: 
the parton-recombination with Q 2 rescaling and the aligned-jet model. Both shadowing 
results are completely different, so that they can be tested by accurate experimental 
measurements. However, both theoretical results are allowed by the existing exper- 
imental data at this stage. The shadowing could be studied at HERA in future by 
observing charged pion productions [8]; hence, there is a good possibility of discrimi- 
nating among various shadowing models by studying the valence-quark distributions. 

4. Nuclear dependence of Q 2 evolution 

The x dependence of nuclear modification in the structure function F 2 has been 
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accurately measured and it is well studied theoretically. On the contrary, there is little 
investigation on Q 2 dependence partly because experimental accuracy was not good 
enough to find nuclear dependence of Q 2 evolution. However, it becomes possible to 
find the details of nuclear Q 2 evolution due to recent NMC analysis of tin and carbon 
F 2 ratios [9]. They found significant deviations d[F$ n /F?]/d[ln Q 2 } ^ 0, so that it is 
an interesting topic for theoretical studies. 

Q 2 dependence of structure functions can be calculated by using the Dokshitzer- 
Gribov-Lipatov-Altarelli-Parisi (DGLAP) equations. In addition, parton-recombination 
(PR) contributions are expected at small x. Although the DGLAP equations are well 
tested by various experimental data, the PR equations are not well established yet. 
There are two possible sources for the nuclear dependence in the evolution equations. 
One is nuclear parton distributions, and the other is the recombination effects. As it 
is discussed in sections 1, 2, and 3, parton distributions are modified due to nuclear 
medium effects. The modifications affect the Q 2 evolution through splitting functions. 
Because parton-recombination probability is proportional to A 1 / 3 , which is the number 
of nucleons in the longitudinal direction, the recombination effects could become sig- 
nificant in large nuclei. We study these two nuclear contributions [10] in comparison 
with the NMC data for d[F^ n / F^]/ d[\nQ 2 }. 

As the input parton distributions in tin and carbon nuclei, we employ those in 
section 1 [1]. They are evolved by leading-order (LO) DGLAP, next-to-leading-order 
(NLO) DGLAP, and PR equations with the help of the computer program bfl.fort77 
in Ref. [11]. Three evolution results are shown at Q 2 =5 GeV 2 together with the pre- 
liminary NMC data [9] in Fig. 5. The DGLAP results agree with the experimental 
tendency; however, the PR curve is well below the data points at small x. The disagree- 
ment between the PR results and the NMC data does not mean the PR mechanism 
is in danger. The large discrepancy is caused first by the choice of initial Q 2 (Qq—0.8 
GeV 2 ) and second by the choice of the constant Kht=1-68 in a higher-dimensional 
gluon distribution. Although the value of Kht is taken according to the Qiu's numeri- 
cal analysis, it is not a uniquely determined quantity. In order to discuss the validity of 
the parton recombination effects, we have to estimate K HT theoretically. From Fig. 5, 
we can at least rule out large recombination contributions. Because the figure indicates 
significant recombination effects, the Q 2 derivative could be used for examining such 
contributions. 

[Summary] The NMC's nuclear Q 2 dependence, d[F$ n / F% ]/<9[ln Q 2 ] ^ 0, could 
be essentially understood by ordinary Q 2 evolution equations together with modified 
nuclear parton distributions. Because the PR evolution results disagree with the data, 
large higher-twist effects from the parton recombinations could be ruled out. However, 
it is encouraging to study the details of the recombination mechanism in comparison 
with the NMC data. 

PS. The Adelaide group for hadron structure (R. P. Bickerstaff, W. Melnitchouk, K. 
Saito, A. W. Schreiber, A. I. Signal, and A. W. Thomas) has made many contributions 
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to similar topics on the nuclear structure functions and on the flavor asymmetry. 



Acknowledgments 

S. K. thanks the Institute of Theoretical Physics at University of Adelaide for its 
financial support for his participating in this conference. This research was partly 
supported by the Grant-in- Aid for Scientific Research from the Japanese Ministry of 
Education, Science, and Culture under the contract number 06640406. 

* Email: kumanos@cc.saga-u.ac.jp. Information on his research is available 
at |http: / / www.cc.saga-u.ac.jp/saga-u/riko/physics/quantuml / structure.html. 



References 

[1] S. Kumano, Phys. Rev. C48 (1993) 2016; C50 (1994) 1247. 

[2] F. E. Close, R. G. Roberts, and G. G. Ross, Nucl. Phys. B296 (1988) 582. 

[3] P. Amaudruz et al. (NMC collaboration), Nucl. Phys. B441 (1995) 3; J. Gomez 

et al. (SLAC-E139 collaboration), Phys. Rev. D49 (1994) 4348. 
[4] S. Kumano, Phys. Lett. B342 (1995) 339. 
[5] S. Kumano, Phys. Lett. B298 (1993) 171. 

[6] R. Kobayashi, S. Kumano, and M. Miyama, Phys. Lett. B354 (1995) 465. 
[7] L. L. Frankfurt, M. I. Strikman, and S. Liuti, Phys. Rev. Lett. 65 (1990) 1725. 
[8] S. Kumano, |nucl-th/ 95 10029] . 



[9] A. Miicklich and A. Sandacz (NMC), private communications on the NMC pre- 

liminary data for d[Fi n / F^]/d[\nQ 2 ]. 
[10] S. Kumano and M. Miyama, |hep-ph/ 95 12244" . 



[11] M. Miyama and S. Kumano, |hep-ph/950824"6| , Comput. Phys. Commun. in 
press; R. Kobayashi, M. Konuma, and S. Kumano, Comput. Phys. Commun. 86 
(1995) 264. 



7 




Fig. 1 x dependence of F 2 /F 2 . 




Fig. 2 Parton recombinations. 
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Fig. 3 Flavor asymmetry in W. 
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